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Abstract 

We study the quantum corrections to the Gross-Pitaevskii equation for two 
weakly linked Bose-Einstein condensates. The goals are: 1) to investigate 
dynamical regimes at the borderline between the classical and quantum be- 
haviour of the bosonic field; 2) to search for new macroscopic quantum coher- 
ence phenomena not observable with other superfluid/super conducting sys- 
tems. Quantum fluctuations renormalize the classical Josephson oscillation 
frequencies. Large amplitude phase oscillations are modulated, exhibiting 
collapses and revivals. We describe a new inter-well oscillation mode, with a 
vanishing (ensemble averaged) mean value of the observables, but with oscil- 
lating mean square fluctuations. Increasing the number of condensate atoms, 
we recover the classical Gross-Pitaevskii (Josephson) dynamics, without in- 
voking the symmetry-breaking of the Gauge invariance. 
PACS: 74.50.+r,03.75.Fi 
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The experimental observation of the Bose-Einstein condensation of a trapped, dilute gas 
of alkali atoms JI], and the high accuracy of the engineering [[2]|§, are opening a new avenue 
to investigate the interplay between macroscopics and quantum coherence. Foundational 
problems of quantum theory and condensed matter |||| can be addressed through 
real (and not just "gedanken") experiments; dynamical regimes not accessible with other 
superconducting/superffuid systems might be testable. 

The main goal of this work is to study the quantum corrections to the classical Gross- 
Pitaevskii dynamics [0 of two weakly linked Bose-Einstein condensates (BEC's) forming a 
Josephson junction. 

The Gross-Pitaevskii equation (GPE) has been shown to describe quite accurately the 
dynamical regimes experimentally investigated so far ||. On the other hand, BEC's can 
be experimentally created with a number of atoms ranging from few thousand to several 
millions, and in a wide variety of confining geometries This is opening the important 

possibility of studying dynamical regimes at the borderline between the quantum and classi- 
cal nature of the bosonic field, and, more generally, to search for new macroscopic quantum 
phenomena not obervables with other superfluid/superconducting systems. 

Different Josephson dynamical regimes are characterized by the ratio of the "Josephson 
coupling energy" Ej and the "on-site energy" E s |9Hl2f . In the limit Ej 3> E s (often referred 



in the literature as "classical" [|T0| , pT| ), both the phase difference and the relative number 
of condensate atoms are well defined. Nevertheless, we will see that quantum corrections 
can significantly modify the classical dynamics even for Ej/E s ~ 10 2 , a regime accessible 
with current BEC technology. Quantum fluctuations renormalize the classical Josephson 
oscillation frequencies. Large amplitude phase oscillations are modulated by partial collapses 
and revivals. It is well known that the relative phase of two weakly linked systems diffuses 
PJT3Hl5| and subsequently revives [[Uj after the suppression of the Josephson coupling. 



In effect, even in two coupled condensates, the relative phase diffuses in the self-trapped, 
running-phase JT^JT^] regime. Then, partial or complete revivals could occur, due to the 
finite number of condensate atoms. An asymmetric potential can also induce phase diffusion 
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In the limit of large condensates, the dynamical equations for the mean values of the 
physical observables decouple from the equations governing the respective quantum fluctu- 
ations, with a smooth crossover from the quantum to the classical GPE regime. 

The classical boson Josephson junction (BJJ) equations, derived by the GPE in the "two- 
mode" approximation |l6|-p0|], can be cast in terms of two canonically conjugate variables: 
the relative population N and phase between the two traps. Quantizing BJJ, the c- 
numbers N and are replaced by the corresponding operators, satisfying the commutation 
relation [0, N] = i |12|^1| . Then the Hamiltonian of two weakly coupled condensates reads 
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H= ^N 2 -Ej cos 4> + AE N (1) 

where Ej (~ a ~ 1) is the "Josephson coupling energy"; E s (~ N^, with (3 = 3/5 
in 3d traps) is the "on-site energy", the analog of the charging energy in a superconduct- 
ing Josephson Junction (SJJ); Nt is the total number of condensate atoms. AE is the 
zero-point energy difference in two asymmetric traps [HJ, or an applied chemical poten- 



tial difference (induced, for example, by the gravitational potential in vertical traps ||). 
The coefficients Ej, E s are determined by the BJJ geometry and the total number of atoms. 
They can be (consistently) calculated as overlap integrals of two orthogonal one-body Gross- 
Pitaevskii wave functions [ |l6|j . 

In the phase representation, the operators in Eq. ([[]) are expressed as N = —2^,0 = 0. 
Then the dynamical equation for the amplitude \I/(0, t) is (h = 1): 



dt 
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r cos ty(4>, t) — iE ^-L, (2) 



where T = En = 2 ^ E ° and the time has been rescaled as — > t. Since we are 
considering an isolated, energy conserving system, the "potential" is periodic and defined 
on a 27r-ring, with the wave-function boundary conditions ^(0) = \]/(0 + 2n). 

In the context of SJJ's, (where only stationary regimes are experimentally accessible), 
Eq. (H) is the drosophila of low-capacitance systems, The effect of dissipation on 



their quantum statistical properties (like phase transitions from normal to superconducting 
phases) has been extensively studied pB| . Other typical effects include the renormaliza- 
tion of the critical Josephson current |24]| , and the macroscopic quantum tunneling among 



metastable minima of the "washboard" potential [PHTTfl. 

In an BJJ, on the other hand, dynamical density oscillations can be studied by shifting 
the position of the laser barrier or tailoring the traps [@,[§. (A similar argument holds when 
considering Raman transitions between two condensate in different hyperfine levels of a 
single traps). The small frequency, < kHz, oscillations of the population imbalance, and the 
corresponding mean square deviations, might be directly monitored by destructive or non- 
destructive techniques. It is worth stressing, then, that the set of experimentally accessible 
observables in BJJ is rather different from the S JJ one. New collective quantum phenomena, 
not accessible with other superconducting/superfluid systems, might be observed in BEC's. 

In this Letter we study analytically the quantum corrections to the classical Gross- 
Pitaevskii- Josephson equations, providing a simple framework to study quantitatively a 
mesoscopic BJJ . 



We consider a time-dependent variational approach. The time evolution of the variational 
parameters is characterized by the minimization of an action with the effective Lagrangian: 

L(ft, ft) =z(W) (3) 

with H = —-§^2 — Tcos0 — iE -^7 and ^(0, ft(t)), ft being the time-dependent parameters. 
This provides the familiar Lagrange equations: 

d dL dL 
dt d(ji dqi 

We choose the time-dependent variational wave function as: 

(0, t) = R[(j) - <j> c (t), X(t)] jStt-MtMQAt)] ( 5 ) 

with the square root of the probability density R and the dynamical phase S being real, 
and: 



S = p(t) [0-0 c (t)] + M [0-0 c (t)] 2 . (6) 

The pairs of time-dependent parameters <f> c (t),p(t) and X(t),5(t) are canonically conjugate: 

OH 

4> c = -^ = 2p + E (7a) 

(7b) 

PjTT 

A=^ = 4« (7c) 
5 = -— = -25 2 + ^-[ RR"dcf>-(V(cf>))] (7d) 

with the effective Hamiltonian: 

H=(T) + (V) = p 2 + 2X5 2 - [* C+n RR"d<\> + (V(<j>)) + £ p. (8) 

J(f>c — ?r 

Thus p(t) is the momentum associated with the center of mass motion <p c (t) = j^-n R 2 {4> — 
4> c ) 4> d(p, and 8(t) is the conjugate momentum of the width of the wave-function, which is 
proportional to A(f) = If^R 2 <p 2 d<f>. The (...) means |^(0 - <f) c )\ 2 ... #, and the 

prime ' stands for J^. The mean value of the population imbalance between the two traps is 
N(t) = {^f(4>,t)N^(4>,t)) = pit), the relative phase is <p(t) = <f) c (t), and the corresponding 
mean square deviations are a%(t) = (N 2 ) — (N) = (T(t)) and cr 2 (t) = 2X(t). 
For r>l, R{4> — <fi c (t), X(t)] can be well approximated by a Gaussian: 

R(M = (4nX)- 1/4 e-^~^ 2 (9) 



with the caveat that during the dynamics its width 2a^,(t) = 2y2X(t) <C 2tt. 
The equations of motion become: 
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iV = -rsin0e-— (10a) 

= 2iV + £ o (10b) 

(7^ = 2(70(5 (10c) 

1 °l 

5 = -25 2 + —j - T cos e~- (lOd) 



with the total (conserved) energy and the relative population dispersion: 



y 4> 
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H = N + a N — T coscj) e~^~ + E N (11) 



aN = 2a~ 4> y 1+A ^ 62 (12) 
The canonically conjugate dynamical variables are N, <fr, as in the classical Josephson Hamil- 



tonian, and the pair -^-,5 — -^Jcr% ~ characterizing the respective quantum fluctua- 
tions. As expected, 0^0$ > | during the dynamics. The classical Josephson equations are 
recovered from Eqs.(|II] a,b) in the limit — > 0. We will discuss more about the transition 
from the quantum to the classical regimes in the following. 
The variational ground state energy of Eq. (D) is given by: 

E 9S = - re~<^ + f = a* , - Te-^ - ^ (13) 

where a^g, ctn,s are the solution of: 

2al e~°l^ = e~^/8al s = I (14) 

The stationary results were first discussed in |Tl|, where Eq.s (p!3|Jl4| ) were obtained min- 



imizing the ground state energy with a Gaussian trial wave function in the case Eq = 0. 
Linearizing Eq. (|10D for small amplitude oscillations, we have: (f> = —2Te~ r7 't''^ 2 (f) The 
condensate atoms oscillate coherently with a frequency (unsealed): 



u g = JE s Ej e~<° /4 (15) 



where the classical Josephson relation gives oj c = ^JE s Ej. The quantum fluctuations renor- 
malize the Josephson plasma frequency, with ^ = exp(— cr^/4) = {a^ s \f2~T)~ l . Notice 
that in the linear regime, the current-phase Eqs.(|T0| a,b) are effectively decoupled from the 
dynamics of the respective fluctuations Eqs.(|T0| c,d). On the contrary, for large amplitude 
4> oscillations, Eqs. ( |I0D cannot be decoupled. In this case the exponential factor modulates 
the amplitude and the frequencies of the oscillations, inducing partial collapses and revivals. 
This can be seen in Fig.(l a-d), where we show the population imbalance, the relative phase 
and the respective mean square deviations as a function of time. 
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Above the critical point (N = 0, <fi = vr/2), the phase <fi(t) starts running, Fig. (2b), and 
the system is set into macroscopic quantum self-trapping (MQST) mode |l6| , |T7|j . The width 



of the wave function grows and the amplitude of oscillations 'collapses', Fig. (2a). In the 
deep MQST regime, when N(t) ~ N(t = 0), the phase diffuses as cr|(i) ~ a\ s + -£i~t 2 , 

4>,s 



Fig(2d), regardless of the initial value of N(t = 0) PTfl . The relative population oscillations 
collapse with a life time r ~ 2Ej 4 E S 4 , while the cr N (t) tends to a constant value, Fig(2c). 
However, since the total number of condensate atoms is finite, the phase can eventually 
revive partially or completely. This can be seen rewriting the wave function in the N 
representation: ^(0, t) = Y,n a N$N(4>) e ~^ ENt > where En are the eigenenergies of Eq. (||). 
For instance, in the limit E s Nt <C Ej, the eigenvalue spectrum is approximately linear 
En ~ , and the revival time is ~ hNx/Ej. More generally, the occurrence of 
a complete or partial revival, or the complete destruction of it, depends on the detailed 
eigenspectrum of the Hamiltonian. We note that Eq. (|1(J) cannot describe the revival after 
a complete collapse since the Gaussian ansatz Eq. @ (and, consequently, the semiclassical 
approximation underlying it) breaks down when ~ 7r. 

Eq.s (|TUD admit, as a dynamical solution, a quite peculiar oscillation mode, with zero 



relative phase and population imbalance, but oscillating fluctuations, according to: 

N(t) = (16a) 
<f>(t) = (16b) 

^) = 4^ + X (16C) 

1 ^ 
(J 4, = — - 2rcr^ e~~ (16d) 

with initial conditions N(t = 0) = 0, <j)(t = 0) = 0; E = and arbitrary a^{t = 0). This 
collective oscillation mode can be experimentally observed by lowering the height of the 
barrier of a BJJ ensemble in thermodynamic equilibrium. This corresponds to changing 
T in Eq. ( |T6| ) from its initial value. The temporal evolution of the ensemble averaged 
observables and the respective mean square deviations can be calculated by tracing the 
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dynamical N(t),<f)(t) trajectories of each junction. 

Classical limit. Increasing the number of atoms, 0$ — > as ~ N T ^ a+l3 \ Moreover, 
for a given initial value N(t = 0), the amplitude of the "particle" oscillations in the <fi- 
potential of Eq. (fj), decreases r-i ~ N T l{a+P) . Then Eqs. (0 a,b) decouple 



from Eqs. ([10] c,d), and the time evolution of the mean values of current and phase become 
independent of the corresponding dispersions. In the MQST regime, the collapse time (and, 
consequently, the time over which the semiclassical predictions are reliable), increases as 
t ~ N^ 313 °^ . In this framework the classical limit emerges naturally, without invoking any 
symmetry breaking argument, as discussed in ||. 

Numerical estimates. Following the analytical estimations of the Josephson coupling 
energy and the on-site energy given for two weakly coupled condensates in [Jl9| , we have: 
T ~ 1 . 7Nt ^ t anh0g/2) ' with a , a s the trap length and the scattering length, respectively, and 
with S ~ j^2ma^(Vo — /i). <Jb is the width of the barrier, Vq its height, and /1 the chemical 
potential. For typical traps and condensates, a ~ 10 4 A, a s ~ 50 A and as ~ 5 jjm. With 
a height of the barrier such that (Vo — /i) ~ 10 nK, we have T ~ 10 — 100 for iVV ~ 1000. 
Varying the width and/or the height of the barrier, and the total number of condensate 
atoms, the system might span from the T 1 to the T ^> 1 limits. The temperature should 
be small compared to the Josephson coupling energy |19| to avoid destroying the quantum 
fluctuations. Damping effects are also reduced by decreasing the total number of atoms. 
Such regimes (vanishing small temperatures, and small population per situ ~ 1000 in an 
optical array) are under current investigations 0. Concluding, we remark that Eq.s (|10"D 
can be straightforwardly generalized to describe interwell tunneling in an array of trapped 
condensates. Work in this direction as well as on the quantitative analysis of the effects of 
temperature and damping on the quantum dynamics are in progress. 

This work was partly supported by NSF Grant PHY94-15583. 
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FIGURES 

FIG. 1. The population imbalance N, the relative phase 4>, and the correspond- 
ing fluctuations and as a function of time. The initial conditions are 
N = 4, (f) = 0, r = 100, 5 = 0, a ( p = 0.26 and E = 

FIG. 2. iV, 0, <Tjv, <T0 as a function of time. The initial conditions are the same as in Fig.(l) 
except N(t = 0) = 50. 
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